(N 

o 

(N 



Origin of magnetic circular dichroism in GaMnAs: giant Zeeman splitting vs. spin 
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We present a unified interpretation of experimentally observed magnetic circular dichroism (MCD) 
in the ferromagnetic semiconductor (Ga,Mn) As, based on theoretical arguments, which demonstrates 
that MCD in this material arises primarily from a difference in the density of spin-up and spin-down 
states in the valence band brought about by the presence of the Mn impurity band, rather than 
being primarily due to the Zeeman splitting of electronic states. 
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Although a multitude of magneto-optical (MO) stud- 
ies have attempted to elucidate the effect of magnetism 
occurring in GaMnAs (both in the epilayer and quantum 
well forms) on its band properties, their results are con- 
flicting [li-ll2|. So far the interpretation of the MO data 
in GaMnAs was carried out using a model similar to that 
used successfully for II-Mn-VI diluted magnetic semicon- 
ductors (DMSs), based on the presence of strong inter- 
actions between magnetic moments and electronic bands 
via s,p-d exchange, that leads to a giant Zeeman split- 
ting of electronic states [13]. Various MO experiments on 
II- VI DMSs yielded a consistent picture of this exchange 
mechanism 13|. In particular, it was found that the p- 
d exchange constant f3 (now known with relatively high 
precision) has antiferromagnetic (AFM) character in all 
II- VI DMSs, due to virtual transitions between the p and 
the half-filled 3d orbitals. 

The same approach, when used to analyze MO data in 
GaMnAs, results not only in a wide range of values for 
/3 @, 0, LL2( but, more importantly, in conflicting reports 
about its sign (FM or AFM) jy, |§, QJ] , raising questions 
about the exchange process between Mn moments and 
band carriers. The correct interpretation of this process 
is of fundamental importance, since it has major conse- 
quences for our understanding of the effect of magnetism 
occurring in any III-Mn-V system on its band proper- 
ties, and therefore for possible future applications. The 
most commonly used technique to study MO properties 
of GaMnAs is magnetic circular dichroism (MCD), which 
arises from magnetic-field-induced differences in the ab- 
sorption of right and left circularly-polarized light, and 
is given by the expression 
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Here T + and T~ are transmission intensities for a + and 
<7~ circular polarizations (which connect the spin-down 
and spin-up holes with the conduction band states, re- 
spectively), and a are the corresponding absorption co- 
efficients. When the difference between a + and a~ near 



the band gap is caused by a Zeeman shift, the polarity 
of the MCD signal identifies the sign of the p-d exchange 
constant. However, the shape and polarity of the MCD 
spectrum in GaMnAs is observed to depend on growth 
conditions [2, lZ|, resulting sometimes in only positive 
MCD signal, and sometimes in both negative and pos- 
itive parts. Application of the Zeeman shift model then 
leads to a FM (positive MCD) or AFM (negative MCD) 
signs for /?, explaining the contradictory claims in the 
literature. 

However, Zeeman shift is not the only mechanism of 
MCD. It can also arise due to magnetic- field-induced 
population changes and/or mixing of states (see, for e.g., 
Ref. Il4 ). In this Letter we present a unified interpreta- 
tion of experimental MCD data, supported by theoretical 
arguments, which demonstrates that MCD in GaMnAs 
arises primarily from a difference in the density of spin- 
up and spin-down states in the valence band, brought 
about by the presence of the Mn impurity band, rather 
than being primarily due to a Zeeman shift. 

GaMnAs samples with thickness of 0.3 /zm were grown 
by molecular beam epitaxy (MBE) on (001) GaAs sub- 
strates. Samples with Mn concentration x ~ 0.01 were 
grown on a 0.4 /xm thick ZnSe buffer, while samples 
with x = 0.0003 and 0.002 were grown on a 0.2 /im 
Gao.70Alo.30As buffer. Both buffers serve as excellent 
etch stops, which allowed us to etch away the GaAs 
substrate in order to eliminate its contribution to op- 
tical transmission. The Mn (and Be, when appropriate) 
concentrations were determined by secondary ion mass 
spectrometry (SIMS), showing a very uniform Mn distri- 
bution throughout the layers. The Curie temperatures 
for samples exhibiting FM behavior were obtained from 
SQUID magnetization data. The MCD measurements 
were performed in the transmission mode using polariza- 
tion modulation via a photo-elastic modulator. 

Figure QJi shows MCD spectra for a Gao.984Mn . i6As 
sample taken in a magnetic field of 5.0 T at tempera- 
tures T = 1.7, 20, 40 and 70 K, i.e. from well below to 
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FIG. 1: a) MCD spectra taken at B = 5.0 T and T = 1.7, 20, 
40 and 70 K on a Gan.984Mno.oi6As layer grown at 280°C. The 
dash-dotted vertical line marks the onset of the sharp rise of 
MCD. b) MCD spectra taken at 5.0 T for similar temperatures 
on a Cdi-zMnzTe sample with x — 0.037. The inset shows a 
schematic illustration of the spin split valence and conduction 
bands in CdMnTe, together with spin conserving transitions 
allowed for a + and o~ circular polarizations (blue and red 
lines, respectively). 



well above its Curie temperature To of 26 K. The MCD 
spectra show a broad positive signal, which rises sharply 
in the vicinity of the energy gap, has a peak around 1.6- 
1.7 eV, and extends well beyond the band gap. Similar 
MCD spectra were observed on FM samples by several 
other authors |l|, |g, [7J. As T increases, the magnitude 
of MCD decreases steadily, while the energy at which 
the sharp rise occurs remains unchanged. This behav- 
ior is drastically different from that observed in II- VI 
DMSs. As an example, Fig. QJ> shows the MCD spectra 
for Cdn.963Mno.o37Te taken at similar values of B and T. 
Here the MCD shows a negative peak (a consequence of 
AFM character of the p-d exchange constant) centered 
at the energy gap. As T increases, the spectrum nar- 
rows as both its onset and its termination energies shift, 
reflecting a strong temperature dependence of the s. p-d 



enhanced giant Zeeman splitting of the band edges [13jj . 
The drastic differences in the shape and temperature de- 
pendence of the MCD spectra of these two materials sug- 
gests that they likely result from different mechanisms. 

The fundamental difference between II-Mn-VI and III- 
Mn-V systems lies in the fact that the divalent Mn acts 
as an acceptor in the III-V host. Compared to Ga 3+ ions, 
Mn 2+ creates a much weaker Coulomb attraction for va- 
lence band (VB) electrons. Taking the host's background 
as reference, this results in an effective repulsive potential 
near the Mn site, pushing one eigenstate above the VB, 
thus forming an impurity state near the Mn ion, while 
the valence band is depleted by one state. Normally this 
state would be empty (i.e., filled by a hole). However, 
in a compensated sample this state can be filled by an 
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FIG. 2: MCD signal (black line) and a~ ,a + (red, blue lines) 
calculated for a ID doped system which is (a) uncompensated, 
(b) partially compensated, and (c) fully compensated. The 
corresponding allowed transitions are shown schematically in 
the insets. The sharp peaks seen at E g reflect the ID density 
of states (DOS) which diverge at the band-edges. 



electron from a donor, such as Asgo ■ Since Mn 2+ is mag- 
netic, the energy of this state is determined not only by 
the effective repulsion, but also by the AFM p-d exchange 
which contributes a fair amount to the total "binding" 
energy [15]. This leads to AFM correlations between the 
spin of Mn 2+ and of the hole (FM correlations with the 
electron, in compensated samples) occupying the impu- 
rity level. Consequently, there is one fewer electron state 
with spin oriented parallel to the Mn spin left in the va- 
lence band. 

We will now discuss the MCD spectrum of a sample 
with a low Mn concentration x, where the impurity lev- 
els have broadened into an impurity band (IB). For that 
purpose we calculated, using the standard Kubo formula 
[161 ] . the a + and a~ absorption spectra using a one- 
dimensional (ID) "toy model". The VB and CB (con- 
duction band) were modeled by nearest neighbor hop- 
ping on two sublattices, with hopping integrals chosen to 
give reasonable effective masses and an on-site potential 
on one sublattice chosen to give a gap of 1.5 eV. Mn 2+ 
ions were assumed to be randomly distributed on sites of 
the VB sublattice. For simplicity, we assumed them to 
be fully spin polarized by an applied magnetic field, and 
chose their impurity potential and AFM p-d exchange 
with VB electrons so as to give a binding energy of 0.1 
eV. The results shown in Fig. [5] are for a system with 
N = 500 sites per sublattice and 10 Mn (a; = 0.02). The 
oscillatory behavior seen at higher energies is a result of 
finite-size effects in the density of states. Minimizing this 
effect requires large N and is the reason why 3D calcu- 
lations, which involve full diagonalization of matrices of 
linear dimension ~ N 3 , are effectively impossible. 

Fig. [5^ shows cv* absorption spectra and the resulting 
MCD signal for a sample with an empty IB (no com- 



pensation). As sketched in the inset, contributions to 
a ± come only from VB — >• CB transitions, both start- 
ing at the band gap energy E g = 1.5 eV. Assuming that 
the Mn spins are fully polarized by the applied magnetic 
field, there are 1 — x times fewer spin-down states than 
spin- up in the VB, and consequently a + absorption is 
weaker than a~ absorption (blue vs. red lines). This 
spin dependent difference in the DOS of the VB leads to 
a positive MCD signal (black line in the figure), which 
rises sharply at the energy gap. If the system is par- 
tially (Fig. [5b) or fully (Fig. [2b) compensated, transi- 
tions from IB to CB will also take place. However, since 
the impurity band states are now occupied by spin-down 
electrons, these transitions will contribute only to the a + 
absorption. Therefore the difference between the two ab- 
sorptions due to IB-to-CB transitions will now result in 
a negative contribution to the MCD signal starting at an 
energy E g — Abe (where Abe is the binding energy), 
i.e., about ~ 1.4 eV, and extending into higher energies. 
Above Eg the MCD signal will have two competing con- 
tributions for compensated samples: a positive contribu- 
tion arising from the difference in the density of states in 
the VB (as discussed above), and a negative contribution 
due to IB-to-CB transitions. The resulting MCD signal 
will then depend on the degree of compensation, as well 
as on the matrix elements of the contributing transitions. 

As discussed above, we used a simplified ID model 
in our calculations. Since the 3D DOS vanishes at the 
band-edges, one expects the MCD signal to peak at 
an energy higher than E g , where the DOS in the VB 
is depleted the most by contributions to the IB. For 
a simple impurity wavefunction (r\I) ~ exp(— r/as), 
the contribution from VB states of momentum k is 
p(k) oc k 2 \(k\I)\ 2 oc k 2 (k 2 + l/a 2 B y 2 (k 2 is due to 
phase-space volume). Contributions from heavy- hole 
(hh) states with energy e^(fc) = h 2 k 2 /2rrih are then 
p(e) = p(k)dk/de ~ ^/e(2mhe/h 2 + \/a 2 B )~ 2 , with a max- 
imum at €h — h /{&rriha B ). Since optical transitions 
conserve momentum, the maximum of the MCD signal 
due to the difference in DOS between the spin-up and 
spin-down states in the VB should then occur at an en- 
ergy E max = E g + e h + e e where e e = ft 2 / (6m* a|) , m* 
being the CB effective mass. For light holes (lh) one 
gets a similar formula, but with the hh mass replaced 
by the lh mass in e/j. Using typical values as ~ 10A, 
ruh ~ 0.5m e , and mi ~ m* ~ 0.07m e gives E max ~ 1.71, 
1.86 eV for hh and lh, respectively. These very rough 
estimates are consistent with the MCD spectrum shown 
in Fig. [Tfej for a sample with very small degree of com- 
pensation. The above model explains not only the sign 
of the MCD signal and the energy range where it peaks, 
but also the temperature independence of the sharp jump 
in the positive contribution to MCD, as seen experimen- 
tally (see Fig. la). According to our model this onset 
depends only on E g , which is temperature independent 
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FIG. 3: Absorption spectra of GaMnAs with x — 0.0003 taken 
at 5 T and T — 2K for circular polarizations <r + and <r~ . The 
corresponding MCD calculated from the a + and a~ absorp- 
tion curves as (a~ — a + )d/2 is also plotted. The MCD is 
dominated by a single broad peak centered at 1.565 eV, cor- 
responding to the energy range where da/dE is negligible. 



in the range studied [17 1 . 

The most direct verification of the proposed scenario 
comes from measurements of a + and aT at B = 5 T 
shown in Fig. [3] for a GaMnAs film with x = 0.0003. 
The sample was grown at a temperature of 600° C in or- 
der to minimize formation of Asca antisites. An antire- 
flection coating was evaporated onto the sample surface 
to suppress Fabry-Perot oscillations for photon energies 
near the band gap. Examination of Fig. 3 indicates that 
the difference between the a + and a~ curves is mainly 
caused by a "vertical" rather than "horizontal" shift, sug- 
gesting that the difference between the two absorptions 
is dominated by the difference in the density of states 
involved in the spin-up and spin-down transitions, rather 
than by the Zeeman shift of the absorption edges. 

We will now compare the predictions of the model with 
the MCD data shown in Fig. [4Jl for several GaMnAs 
films. Samples A, B and D have the same Mn con- 
centration x — 0.014, but were grown under different 
conditions, resulting in different levels of compensation: 
Sample D was grown at 200°C; samples A and B were 
grown at 250°C; and sample A was additionally co-doped 
with Be. As a result, sample D is expected to be highly 
compensated, while sample A has the lowest degree of 
compensation. As our model predicts, the below-gap 
negative contribution to MCD signal is the strongest in 
sample D, still visible in sample B, but absent in sam- 
ple A. Above the band gap the positive MCD signal for 
the highly-compensated Sample D is considerably weaker 
than for samples A and B. This again is consistent with 
our model, because the above-£" g MCD signal has two 
opposite contributions; and since in Sample D the num- 




FIG. 4: a) MCD spectra taken at T = 1.8 K and B = 5.0 
T for several GaMnAs samples. Samples A, B, and D have 
x = 0.014, while sample C has x = 0.002. For other details, 
see text, b) MCD spectra taken on sample B at T — 1.8 K 
and B = 0, 0.5, 2 and 5 T. The horizontal line marks MCD 
= 0. Inset: Magnetic field dependence of MCD for sample B 
observed at T = 1.8 K at photon energy 1.405 eV, showing a 
small but finite hysteresis loop. 



ber of spin-down electrons occupying IB is roughly equal 
to the number of spin- up states depleted from the VB, 
the resulting signal should be significantly smaller than 
in samples with lower compensation. 

We now consider the shape of the above-gap spectrum. 
Sample B provides clear evidence for the existence of two 
MCD peaks above E g , at locations consistent with our 
rough estimates for hh and lh contributions. This identi- 
fication is further supported by the presence of only one 
peak (at ~1.65 eV) in sample C, grown at T$ — 260°C 
with a much lower Mn concentration x — 0.002. Since 
in the single-impurity limit the largest contribution to 

3 



we 



the impurity state comes from heavy holes 
expect to see only the hh peak in the MCD signal for 
very low x, as indeed evidenced by sample C. However, 
as x increases and the hh states around the optimum en- 
ergy become exhausted, one expects to see more of the 
lh states as well as states deeper in the VB being pulled 
into the IB. This is consistent with the emergence of the 
second peak (around 1.85 eV) in sample B, and the broad 
unresolved peak at ~1.8 eV in sample A. 

A final check of our scenario comes from the variation 
of MCD with changing magnetic field, shown in Fig. 0t> 
for sample B. The FM character of this sample (Tq = 16 
K) is evidenced by the presence of MCD signal at zero 
magnetic field, where, however, we only observe the pos- 
itive part of MCD produced by interband transitions. 
The absence of the negative contribution at B = indi- 
cates that, even though the sample is FM, the IB states 



filled by electrons are not spin-polarized at B = 0. Only 
the empty (i.e. hole-occupied) states in the IB are spin- 
polarized, thereby giving rise to a difference in the spin- 
up and spin-down DOS in the VB, and consequently to 
a positive MCD signal above E g . This is fully consistent 
with the picture proposing that the Mn-Mn interaction in 
low-x samples is mediated by holes hopping inside an IB 



21]. Only the Mn ions with high overlap with these 
holes can be ordered magnetically at Tc, and in turn 
polarize the spins of the holes weakly bound to them. 
The Mn spins with low overlap with the holes, i.e., those 
with impurity states filled with electrons, remain para- 
magnetic down to very low temperatures. Since these 
Mn are not polarized, neither are the electrons filling the 
states bound to them. 

When a magnetic field is applied, however, it polarizes 
these Mn spins and they in turn polarize their corre- 
sponding bound states, which results in a gradual emer- 
gence of IB-to-CB transitions, and thus in an onset of 
a negative contribution to MCD below E g . This is pre- 
cisely the behavior seen experimentally in Fig. @J). What 
follows then is that the negative part of MCD must mir- 
ror the paramagnetic behavior of the relevant Mn spins. 
This is verified in the inset in Fig. [4]b, which shows that 
the magnetic field dependence of the MCD signal mea- 
sured at 1.4 eV follows a paramagnetic behavior, even 
though the hysteresis loop seen in the inset indicates a 
FM character of the sample as a whole. 

In conclusion, this work explains in a unified manner 
the complex behavior of MCD in GaMnAs that has long 
been a puzzle in this field. The solution is drastically dif- 
ferent from that based on the giant Zeeman shift relevant 
to II-Mn-VI DMS, and underscores again the importance 
of the impurity band to our understanding of III-Mn-V 
DMSs. 
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